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Various regularization methods have been used to compute the self-force acting on a static particle 
in a static, curved spacetime. Many of these are based on Hadamard's two-point function in three 
dimensions. On the other hand, the regularization method that enjoys the best justification is that 
of Detweiler and Whiting, which is based on a four-dimensional Green's function. We establish the 
connection between these methods and find that they are all equivalent, in the sense that they all 
lead to the same static self-force. For general static spacetimes, we compute local expansions of the 
Green's functions on which the various regularization methods are based. We find that these agree 
up to a certain high order, and conjecture that they might be equal to all orders. We show that this 
equivalence is exact in the case of ultrastatic spacetimes. Finally, our computations are exploited to 
provide regularization parameters for a static particle in a general static and spherically-symmetric 
spacetime. 



I. INTRODUCTION 

A test body moving freely in a curved spacetime follows 
a geodesic of the spacetime. When, however, the body 
carries a (scalar or electric) charge, the field created by 
the charge interacts with the spacetime curvature in such 
a way as to produce a deformation of the field lines from 
an otherwise isotropic distribution around the body. The 
field gives rise to a net self-force acting on the body, and 
the self-force prevents it from moving on a geodesic. The 
self-force typically contains two components, a radiation- 
reaction force that is accompanied by a loss of energy 
to radiation, and a conservative force that survives even 
when the body is maintained in a stationary position. A 
self-force can also be present in the absence of a charge, 
when the body's mass is too large for it to be considered 
a test mass; in this case the body creates a gravitational 
perturbation that affect its motion, which is no longer 
geodesic in the background spacetime. The (scalar, elec- 
tromagnetic, and gravitational) self-force has been the 
topic of intense development in the last several years; for 
an extensive review see Ref. [l[ . Most of this activity was 
focused on the gravitational case, in an effort to model 
the inspiral and gravitational- wave emissions of a binary 
system with a small mass ratio 0-0]. 

Self-force computations are usually attempted under 
the assumption that the body is a point particle, in order 
to avoid the largely irrelevant complications associated 
with internal structure. In this context, however, the 
very definition of the self-force requires scrutiny. Given 
that the field of a point particle diverges at the position 
occupied by the particle, it is not immediately clear how 
one can make sense of its action on the particle and con- 
struct a self-force that is well defined, finite, and in agree- 
ment with the self-force acting on an extended body in 



the the limit in which the size is taken to zero. One must 
find a sensible regularization procedure that not only re- 
turns a finite expression for the self-force, but does so in a 
unique and physically well-motivated way. In this paper 
we examine regularization procedures that have been in- 
voked in the computation of (scalar and electromagnetic) 
self-forces in the restricted context of static particles in 
static spacetimes. Our aim is to show that the differing 
procedures are equivalent and lead to the same self-force. 
To the best of our knowledge, this issue has not been pre- 
viously addressed in the literature. 

In our view, the regularization procedure that has re- 
ceived the best physical and mathematical justification 
is the one proposed by Detweiler and Whiting [f|. The 
method, which is completely general and not restricted 
to static situations, involves a decomposition of the field 
created by the particle into singular and regular pieces. 
The singular field is precisely identified by a local con- 
struction, and is designed to provide an exact solution 
to the field equation sourced by the particle, with the 
property that it shares the singularity structure of the 
particle's actual field. The regular field is the difference 
between the actual field and the singular field; it satisfies 
the source-free field equation, it is smooth at the parti- 
cle's position, and it is known to be entirely responsible 
for the self-force. The Detweiler- Whiting regularization 
method has been thoroughly justified 0-0]' an d ^ nas 
emerged as the method of choice in most self-force com- 
putations reported in the recent literature. Because of its 
generality and naturalness, it is the standard by which 
other regularization methods must be compared. 

Many self-force computations, however, did not make 
use of the Detweiler- Whiting regularization procedure, 
but employed instead ad hoc procedures that perhaps do 
not enjoy the same degree of justification. This is the case 
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of all computations of self-forces acting on static parti- 
cles in static spacetimes [Tol - [2"4| , which involved a variety 
of regularization methods. In the pioneering Smith- Will 
paper [lfj, for example, the field of a static electric charge 
in the spacetime of a Schwarzschild black hole was reg- 
ularized by the Copson solution (25|, which was shown 
to be as singular as the particle's own field but to exert 
no force. As other examples, self-force computations for 
charges in wormhole s pac etimes [l7l - l20| . or for charges 
near global monopoles [2l| - [24j , were regularized with the 
help of Hadamard's two-point function, defined in each 
spatial section of the four-dimensional spacetime (or in 
a conformally related space). Because Copson's solution 
is known to be an exact representation of Hadamard's 
function in the (conformally related) spatial sections of 
the Schwarzschild spacetime, these regularization meth- 
ods are essentially the same. 

The issue that interests us in this paper is the re- 
lationship between these regularization procedures, and 
whether they can be shown to be equivalent, so that they 
will lead to the same self-force. The regularization pro- 
cedures mentioned previously are all based on a choice of 
Green's function for the (scalar or electromagnetic) field. 
We shall consider a number of possible choices. 

The first is the four-dimensional version of the 
Detweiler- Whiting singular Green's function, given by 

G\{x, x') = ±U(x, x')S(a) - ±V(x, x')e(a), (1.1) 

in which x and x' are spacetime events, assumed to be 
sufficiently close that they are within each other's nor- 
mal convex neighborhood, a := a(x,x') is Synge's world 
function, equal to half the squared geodetic distance be- 
tween x and x' , O is the Heaviside step function, 5 is the 
Dirac distribution, and U, V are two-point functions that 
are known to be smooth when x — > x' . Because a = 
when x and x' are linked by a null geodesic, we see that 
the Green's function is singular on the past and future 
light cones emerging from x' , and has support outside 
the light cones, where a > 0; it is also symmetric in its 
arguments. As stated previously, the Detweiler- Whiting 
Green's function gives rise to a robust regularization pro- 
cedure that applies to any particle moving in any space- 
time. 

For static particles in static spacetimes, an adequate 
substitute for the four-dimensional singular Green's func- 
tion is its three-dimensional variant 



by 



G s 3 (x,x') 



Gl(x,x')dr', 



(1.2) 



obtained by integrating the four-dimensional Green's 
function over the proper time t' of a static observer at the 
spatial position x'. For static particles in static space- 
times, this Green's function gives rise to the same regu- 
larization procedure as the four-dimensional version. 

An alternative choice of Green's function, which is also 
appropriate in the case of static particles in static space- 
times, is the three-dimensional Hadamard function, given 



(1.3) 



in which <r(x, x') is now half the squared geodetic dis- 
tance between x and x' as measured in the purely spa- 
tial sections of the spacetime, and W is a smooth two- 
point function. This Green's function, also known as 
Hadamard's elementary solution (2(| , is a local construc- 
tion that is known to reproduce the singular behavior of 
a field sourccd by a point particle at x' . It is a plausi- 
ble starting point for a regularization procedure, but as 
stated above, it does not enjoy the same level of justifi- 
cation as the Detweiler- Whiting Green's function. 

As a final choice we shall also consider a Green's func- 
tion Gg(a;,a;') that is related to the Hadamard function 
by a conformal transformation. This is to account for 
the fact that it is often convenient, when solving for the 
field in the spatial sections of the spacetime, to formulate 
the field equation in a conformally related space. If the 
metric on the original spatial sections is h a b, then the 
metric on the conformally related space is h a t, = Q~ 2 h a b, 
in which £l(x) is a scalar field. When h a b is simple, the 
field equation simplifies in the conformally related space, 
and the field can then be regularized with the help of G3 , 
which differs from G3 by factors of f2. 

Our main goal in this paper is to compare the reg- 
ularization procedures that are based on Gf(x,x'), 
G^(x,x'), and G^(x,x'). Our preferred regulariza- 
tion method is the one based on the Detweiler- Whiting 
Green's function, because the resulting singular field was 
proved to share the same singularity structure as the par- 
ticle's actual field and to exert no force on the particle. In 
this case we find that if e is a measure of distance between 
x and x', then the singular Green's function admits the 
local expansion 



G s 3 (x,x f ) 



l + .9ie + .9 2 e 2 +.9 3 e 3 + 0(e 4 ) , (1.4) 



with expansion coefficients 171, g2, and 33 (computed be- 
low) that depend on geometrical quantities (such as the 
spatial Ricmann tensor) evaluated at x' . Removing this 
from the particle's actual field returns a regularized field 
that is finite (indeed, once diffcrentiable) at the position 
of the particle, leading to a straightforward computation 
of the self-force. 

Our main result is the statement that the local expan- 
sions of G^(x,x') and G^(x,x') are identical to the lo- 
cal expansion of the Detweiler- Whiting Green's function, 
and that they therefore lead to equivalent regularization 
procedures. We view this as a significant result that clar- 
ifies and justifies the alternative regularization methods 
that have been employed in self-force computations; it 
puts these computations on a firmer footing, and lends 
them additional credence. 

We show by explicit calculation that the local expan- 
sions of all Green's functions agree through order e 2 , 
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which is more than enough to guarantee the same self- 
force, but we also ask whether equality could be estab- 
lished to all orders in e. We provide only a partial answer 
to this question. We collect evidence that the conjectured 
equality is likely to be true, sketch a proof that relies 
on a strong assumption about the convergence of formal 
power series, and present a complete proof of equality in 
the special case of ultrastatic spacetimes. 

As an application of our results we consider the scalar 
and electromagnetic self-force acting on a particle held 
in place in a static and spherically-symmetric space- 
time. The assumption of spherical symmetry implies that 
the self-force can be easily computed with a mode-sum 
method based on a spherical-harmonic decomposition of 
the field. The mode coefficients of the singular field then 
play the role of regularization parameters [27T - l29l | that 
can be inserted in the mode-sum to ensure convergence. 
We use our local expansions to compute these regulariza- 
tion parameters for any static and spherically-symmetric 
spacetime. We expect that our explicit listing of regular- 
ization parameters will greatly facilitate future self-force 
computations. 

Our developments below rely heavily on the general 
theory of bitensors and Green's functions in curved space- 
time, as developed in Ref. [3(| and summarized in Ref. [![, 
from which we import our notations. We begin in Sec. |TT] 
with a description of the static spacetimes that are impli- 
cated in this work. We continue in Sec. IIIII with a review 
of the scalar and electromagnetic field equations in static 
spacetimes, along with the associated Green's functions, 
and in Sec. IIVI we review the Hadamard construction of 
the three-dimensional Green's functions. The following 
sections contain our new work. In Sees. |V] and I VII we 
compute the local expansions of the Hadamard functions. 
We do the same for the Detweiler- Whiting functions in 
Sec. lVIll and prove the equality of the local expansions in 
Sec. I Villi In Sec. IIXI we prove that in ultrastatic space- 
times, the Hadamard and Detweiler- Whiting functions 
are strictly equal to one another. And finally, in Sec. |X] 
we consider the scalar and electromagnetic self-forces in 
static, spherically-symmetric spacetimes, and involve the 
local expansions in a computation of regularization pa- 
rameters for mode-sum computations of the self-force. 

In the following we denote a spacetime event by x or 
x' , and a spatial position by x or x' , so that x = (t,x) 
and x' = (t',x'). Spacetime tensors at x are denoted 
A a , with a Greek index a that ranges over the values 
{0, 1, 2, 3}; spacetime tensors at x' are denoted A a , with 
a primed Greek index. Spatial tensors at x are denoted 
A a , with a Latin index a that ranges over {1, 2, 3}; spatial 
tensors at x' are denoted A a , with a primed Latin index. 
Variants of these notations will be introduced as needed. 



II. STATIC SPACETIMES 

The class of spacetimes considered in this paper admits 
a hypcrsurface-orthogonal, timelikc Killing vector t a , and 



in an adapted coordinate system the metric is expressed 
as 

ds 2 = -N 2 dt 2 + h ab dx a dx b , (2.1) 

in terms of a lapse function N and a spatial metric h ab 
which depend on the spatial coordinates x a only. No 
other assumptions are placed on the spacetime. We in- 
troduce the vector field 

A a :=d a \aN=^, (2.2) 

which acts as a substitute for d a N; the vector has a van- 
ishing time component. 

A straightforward computation reveals that the con- 
nection coefficients are given by 

4 Tl = A ai *T? t =N 2 A a , 4 T a bc = T a bc , (2.3) 

in which T^ c is the connection compatible with the spatial 
metric h ab . We shall indicate covariant differentiation rel- 
ative to the spacetime connection with the operator V p 
or with a semicolon (for example, V^A" — A"), and 
covariant differentiation relative to the spatial connec- 
tion with the operator D a or with a vertical stroke (for 
example, D a A b = A b .). 

For future reference we examine a vector v a that 
is parallel-transported along a purely spatial curve de- 
scribed by the parametric equations t — constant, x a = 
z a (s), in which s is proper distance. The vector tan- 
gent to the curve is n a with components n* = and 
n a = dz a /ds. The equation of parallel transport is 

di) a 

^ + 4 r^/V = o, (2.4) 

and its time component reduces to dv* / 'ds+(A a n a )v t = 0. 
With A a n a = N~ l dN/ds, the solution to the differential 
equation is Nv l = constant, or 

vt{s) =W) vt{Q) - (2 - 5) 

For the spatial components we find that the equation 
reduces to dv a /ds + T b l c v b n c = 0, which states that v a 
is parallel-transported as if it were a vector in a three- 
dimensional space with metric h ab . When the spatial 
curve is a spacetime geodesic, we find that Eq. (|2.5[) 
produces n'(s) = 0, while the spatial components re- 
veal that n a satisfies the geodesic equation in the three- 
dimensional space. A purely spatial geodesic in space- 
time is therefore a geodesic in a three-dimensional space 
with metric h ab . 

The nonvanishing components of the Ricmann tensor 
are 

A Rtatb = N 2 (A a \ b + A a A b ^j, A Rabcd = Rabcd, (2-6) 

in which R a bcd is the Ricmann tensor associated with the 
spatial metric h ab . The nonvanishing components of the 
Ricci tensor are 

4 R tt =N 2 (A c lc + A c A c ) 7 (2.7a) 
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A Rab = Rab - (4.|6 + 44) ■ (2.7b) 
The Ricci scalar is 

4; 



l R = R-2(A c ]c + A c A 



(2.8) 



For future reference we also record the components of the 
covariant derivative of the Riemann tensor: 

A Rtatb:c = N 2 (4, bc + 4, c 4 + 44|c) , (2.9a) 

A Rtabc-t = N 2 (4|&4 ~ 4|c4 - A d R dabc ) , (2.9b) 



Rabcd;e — Rabcd\e 

We also have 



(2.9c) 



R tt , a = N 2 (A c [ca + 2A c A c]a ) , (2.10a) 
Rta-t = -N 2 (A c {c A a - A c {a A c + A c R ca ) , (2.10b) 



A 



a\bc 



4|c4 — 44|c 



and 



l R ia = R la -2(A c lca + 2A c A cla 



(2.10c) 



(2-11) 



for the covariant derivative of the Ricci tensor and scalar. 

Below we shall consider a scalar or electric charge at 
rest in the static spacetime. The charge follows an orbit 
of the timelike Killing vector, and the only nonvanishing 
component of its velocity vector is 



1 

N' 



(2.12) 



The covariant acceleration is defined by a" := u a .^u^, 
and its nonvanishing components are 



a a = A a . 



(2-13) 



We shall also require the vectors a a := a a .pvr and a" := 
h a .Q%iP\ the nonvanishing components are 



6 1 = j^A c A c , a a = (A c A c )A a . 



(2.14) 



III. FIELD EQUATIONS AND GREEN'S 
FUNCTIONS 



A. Scalar field 

The potential $ generated by a scalar-charge density 
fi obeys the wave equation 



(3.1) 



in any four-dimensional spacetime; □ := <? Q/3 Vo,V^ is the 
covariant wave operator. Our considerations in this pa- 
per are limited to scalar fields that are minimally coupled 
to the spacetime curvature; it is, however, a very straight- 
forward exercise to extend our discussion to arbitrary 



couplings. When the spacetime is static, and when the 
potential and charge density are both time- independent, 
the wave equation reduces to 



V 2 $ + A a d a <i> = -Ann; 



(3.2) 



here V 2 := h ab D a D b is the covariant Laplacian operator 
in a three-dimensional space with metric hab- 
it is sometimes convenient to formulate Eq. (|3.2[) in a 
conformally related space with metric h ao ; this is related 
to the original metric h a b by the conformal transforma- 
tion 



fl 2 h 



ab i 



(3.3) 



in which Q is a function of the spatial coordinates x a . 
As a consequence of this transformation we find that 
h ab = n- 2 h ab and h 1 ' 2 = n 3 h^ 2 . A simple computa- 
tion reveals that in the conformal formulation, Eq. (|3.2[) 
becomes 



V 2 $ + i a d Q $ = -47T/i, 



(3.4) 



in which V 2 := h ab D a Db is the covariant Laplacian op- 
erator associated with the metric h a b, 



A a := h ab d b ]n(Nfl) = h 



a b 



d b N 
N 



and 



ft := il 2 /i. 



(3.5) 



(3-6) 



The four-dimensional Green's function associated with 
Eq. (|3.ip is G^(x, x'), which satisfies the wave equation 



\3Gi(x, x) = — 47r<54(x, x), 



(3.7) 



in which S4,(x, x') is a scalarized Dirac distribution de- 
fined by 



di(x, x') 



S(x-x') 



(3.8) 



where S(x — x') is the usual product of four coordinate 
delta functions. A solution to Eq. (|3.1j) is then given by 



= / G 4 (x,x')fj,(x') \T^g'd 



4 / 



(3-9) 



The three-dimensional Green's function associated 
with Eq. (|3.2[) is G3(x,x'), which satisfies the Poisson 
equation 

V 2 G 3 (a;, x') + A a d a G 3 (x, x') = -4ttS 3 (x, x'), (3.10) 

in which 5^{x,x') is a scalarized Dirac distribution de- 
fined by 



5 3 (x,x') 



6(x-x') 



(3.11) 
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where 5(x — x') is the usual product of three coordinate 
delta functions. A solution to Eq. (|3.2p is then given by 



G 3 (a;,3;')/i(a;')VW 3 3;' 



(3-12) 



The Green's function associated with Eq. p.4[) is 
G 3 (x,x'), which satisfies 

W 2 G 3 (x, x') + A a d a G 3 (x, x') = -4ttS 3 (x, x'), (3.13) 

in which 5 3 (x, x') is defined by 

S(x - x') 



6 3 (x,x') 

A solution to Eq. (|3.4p is then given by 

$(x)= / G 3 (x,x')p,(x') Vh/d 3 x'. 



(3.14) 



(3.15) 



The relation between G4 and G 3 can be identified by 
performing the time integration in Eq. (|3.9[) and compar- 
ing with Eq. (|3~T2|) . The result is 



G 3 (x,x')= / G 4 (x,x')N(x')dt'. 



(3.16) 



Apart from the factor of N(x'), which converts from co- 
ordinate time to proper time at x' , the three-dimensional 
Green's function is simply the time integral of the four- 
dimensional Green's function. 

The relation between G 3 and G 3 is found by making 
the substitutions h 1 / 2 = Q^h 1 / 2 and ji = fl~ 2 fi within 
Eq. (|3.12[) and comparing with Eq. (|3.15[) . The result is 



G 3 (x,x') = n(x')G 3 (x,x'). 



(3.17) 



This can be confirmed by expressing Eq. (|3.10[) in terms 
of the conformally-related metric h a t,. We find that the 
equation becomes 



V 2 G 3 (x, x') + A a d a G 3 (x, x') = -4tt 



5{x - x') 



QWh 



, (3-18) 



and comparison with Eq. (|3.13[) allows us to make the 
identification of Eq. (|3.17p . 

When the charge density /x describes a point charge 
q moving on a world line 7 described by the parametric 
relations z(t), we have that 

fx(x) = q / 6 4 (x,z(t)) dr. (3.19) 
J j 

For a general world line the scalar charge produces a 
potential given by Eq. (|3.9[) . which evaluates to 



$(x) = q / G 4 (x, z(t)) dr. 
'1 



(3.20) 



For a static charge at a fixed position z, the integral of 
Eq. (|3.19[) evaluates to 



H(x) = q6 3 (x,z), 



(3.21) 



and in this case the potential, as given by Eq. (|3.12D . 

becomes 



$(sc) = qG 3 (x,z). 



(3.22) 



The link between Eqs. (|3.20j) and p.22p can be seen di- 
rectly from Eq. ()3.16[) . 

In the conformal formulation we have instead 



fi(x) = q 



S 3 (x,z) 

n(z) ' 



and Eq. (j3.15[) produces 



= q 



G 3 (x,z) 

n(z) ' 



(3.23) 



(3.24) 



This result is compatible with Eq. (|3.22p by virtue of 
Eq. (13171) . 



B. Electromagnetic field 

A current density j a creates an electromagnetic field 
F a p that satisfies Maxwell's equations 



F° 



[p = Awj a , F aPn + F ia , p + F M , a = 0. (3.25) 



The homogeneous equations are automatically satisfied 
when the field tensor is expressed in terms of a vector 
potential $ Q , 



F afj = V a <f> f3 - V/3$ c 



(3.26) 



The inhomogencous equations then take the form of a 
wave equation for the vector potential, 



(3.27) 



provided that $ a is required to satisfy the Lorenz gauge 
condition 



0. 



(3.28) 



In a static spacetime, and for a static distribution of 
charge, the only relevant component of Maxwell's equa- 
Anj*, and with F ta := —d a $>t this ra- 



tions is F 



duces to 

V 2 $ 4 - A a d a <S> t = 4tt M , h := N 2 f = -j t . (3.29) 



This equation also follows from evaluating Eq. (|3.27p in 
a static spacetime. Comparing with Eq. (|3.2p , we see 
that (f> t satisfies the same Poisson equation as a scalar 
potential except that the sign of A a is reversed. It is 
easy to see that the gauge condition of Eq. (|3.28p becomes 
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D a <& a + A a $ a = in a static spacetime, and that it has 
no impact on $ 4 . 

In a conformal formulation in which the spatial metric 
is expressed as h a b = f2 2 /i a 6, Eq. (j3.29j) becomes 

V 2 $ t - A a d a ®t = 47r/i, (3.30) 

in which 

(3 . 31) 

and 

p, = fl 2 pi = N 2 Q 2 f . (3.32) 

The four-dimensional Green's function associated with 
Eq. (|3.27[) is G£ (x, x'), which satisfies the wave equation 

nGf(x,x')-R^Gf(x,x') = -iTrgfix^y^x'), 

(3.33) 

in which g£ (x,x') is an operator of parallel transport, 
taking a vector Api at x' and producing a parallel- 
transported vector A a at x. A solution to Eq. (|3.27[) 
is then given by 

* a (s)= [G/(x,xr)jptf)y/=g'd*xr. (3.34) 

The three-dimensional Green's function associated 
with Eq. (|3.29|) is G 3 (x,x'), which satisfies the Poisson 
equation 

V 2 G 3 {x, x') - A a d a G 3 (x, x') = -iirS a {x, x'). (3.35) 
A solution to Eq. (|3.29j) is then given by 

$ t (aj) = - J G 3 (x, x')n{x') Vh/d 3 x'. (3.36) 

Notice the minus sign on the right-hand side of Eq. (|3.36|) , 
which can be compared with its scalar equivalent in 
Eq. (|3.12[) . Notice also that while we denote both 
the scalar and electromagnetic Green's functions by 
G 3 (x,x'), these functions are not equal to each other 
because they satisfy distinct differential equations. 

The Green's function associated with Eq. p.30|) is 
G 3 (x, x'), which satisfies 

V 2 G 3 (cc, x') - A a d a G 3 (x, x') = -Att~5 3 {x, x'). (3.37) 

A solution to Eq. (|3.30|) is then given by 

$ t (sc) = - J G 3 (x, x')p,(x') Vh/d 3 x', (3.38) 

which features the same minus sign as in Eq. p.36|) . 

The relation between G£ (x,x') and G 3 (x 7 x') can 
be identified by performing the time integration in 
Eq. (|3.34p and noticing that in a static situation, the 



integral involves j t > only. Comparing with Eq. (|3.36[) pro- 
duces 

G 3 (x,x') = J G t t '(x,x r )N(x')dt', (3.39) 

essentially the same relation as in the scalar case. 

The relation between G 3 and G 3 is found by making 
the substitutions h 1 / 2 = Vfih 1 / 2 and /j = Q~ 2 jl within 
Eq. (|3.36p and comparing with Eq. (|3.38[) . The result is 

G 3 (x,x f ) = n(x')G 3 (x,x'), (3.40) 

the same relation as in the scalar case. 

When the current density j a describes a point charge 
e moving on a world line 7 described by the parametric 
relations z(r), we have that 

j a {x) = e [ g^(x, z)u»6 4 (x, z) dr. (3.41) 

For a general world line the electric charge produces a 
potential given by Eq. (|3.34p . which evaluates to 

$ a (x) = e I G a ^{x,z)u^ dr. (3.42) 

For a static charge at a fixed position z, the integral of 
Eq. (|3.4ip evaluates to f{x) = eN^ 1 (z)S 3 (x,z), so that 

fi{x)=eN(z)5 3 {x,z). (3.43) 

In this case the potential, as given by Eq. (|3.36p . becomes 

$ t (x) = -eN(z)G 3 (x,z). (3.44) 

Notice the extra minus sign and factor of N when com- 
paring this with Eq. (|3.22p . In the conformal formulation 
we have instead 

m = e^-~8 3 {x,z), (3.45) 
and Eq. p.38p produces 

N(z) ~ 

^t(x) = -e-^G 3 {x,z). (3.46) 
This result is compatible with Eq. (|3.44|) by virtue of 

Eq. ngm 

IV. HADAMARD'S CONSTRUCTION 
A. Scalar field 

We wish to find a representation for a Green's function 
G 3 (x,x') that satisfies 

W 2 G 3 (x,x') + A a d a G 3 (x,x') = -4n6 3 {x,x'). (4.1) 
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The Hadamard's construction developed here applies to 
this equation — a copy of Eq. (|3.10p — but it applies 
just as well to the conformally-related formulation of 
Eq. (13~T3l . 

The Hadamard construction for the Green's function 

is m 

<#(*,*0 = -^£> (4.2) 

in which the two-point function W(x,x') satisfies the 
differential equation 

2a a d a W+(V 2 a+A a a a -3)W~(2a)(V 2 W+A a d a W) = 

(4.3) 

together with the boundary condition 

W(x',x') = l: (4.4) 

here a a := da/dx a and V 2 er := h ah D a <Jb. The two-point 
function is known to be smooth in the coincidence limit 
x — > x', so that the factor (2ct) -1 / 2 is fully responsi- 
ble for the singular behavior of the Green's function at 
coincidence. 

To construct W we express it as an expansion in powers 
of 2<j, 

OO 

W(x,x') =J2w n (x,x f )[2a(x,x f )] n , (4.5) 

insert this within Eq. (|4.3[) . and collect powers of 2a, 
making use of the identity a a a a = 2a. Setting each co- 
efficient to zero, we find that each W n must satisfy the 
differential equation 

V 2 W^_! + A^aWn^ = 2(1 - 2n)a a d a W n 
+ (l-2n)(V 2 a + A a a a )W n 
- [3 + 4n(n-2)]W„. (4.6) 

The burden of enforcing Eq. (|4.4p is then placed solely 
upon Wo, which must satisfy 

W Q (x',x') = l. (4.7) 

Equation (|4.6[) is a recursion relation for each W n . With 
W n -i previously determined, W n is obtained by select- 
ing a base point x' and integrating Eq. (|4.6[) along each 
geodesic that emanates from x'. 

B. Electromagnetic field 

We now wish to find the Hadamard representation for 
the electromagnetic Green's function G 3 (x,x'), which 
satisfies 

W 2 G 3 (x,x') - A a d a G 3 (x,x') = -4ir6 3 (x,x'). (4.8) 

The Hadamard construction applies to this equation — 
a copy of Eq. f|3 . 35[) — but it applies just as well to the 
conformally-related formulation of Eq. (|3,37[) . 



The construction is obtained directly from the scalar 
case by altering the sign of A a in all equations. The 
Hadamard representation for the Green's function is 

G»{x,x>) = ^^ (4.9) 
y/2a(x, x') 

in which the two-point function W[x, x') is smooth in 
the coincidence limit x — >• x' . It satisfies the differential 
equation 

2a a d a W+(V 2 a-A a a a -3)W-(2a)(V 2 W-A a d a W) = 

(4.10) 

together with the boundary condition 

W(x',x') = l. (4.11) 

As in the scalar case we express W as an expansion in 
powers of 2a, 

OO 

W(x, x') = ]T W n (x, x') [2a{x, x')] (4.12) 

in which each coefficient W n must satisfy the differential 
equation 

V 2 W„_i - A a d a W n - X = 2(1 - 2n)a a d a W n 
+ (l-2n)(V 2 a-A a a a )W n 
- [3 + 4»(n-2)]W„. (4.13) 

V. LOCAL EXPANSION OF HADAMARD'S 
FUNCTION 

A. Scalar field 

We wish to express the three-dimensional Green's func- 
tion G3 (x, x') as a local expansion about the base point 
x' . We return to the Hadamard construction of Eq. (|4.2[) 
with the expansion of Eq. (|4.5p , and now express Wq and 
W\ as the local expansions 

Wo = 1 + W°a a ' + \w%a a ' ' a h ' + \w° a , Vc ,a a ' ' a v >' 
+ 0(e 4 ) (5.1) 

and 

Wi = W 1 + W x a ,a a ' + 0(e 2 ), (5.2) 

in which a a t := da/dx a and each expansion coefficient 
is an ordinary tensor at x'. We let e be a measure of 
the distance between x and x' , and the expansions of 
Eqs. (|5.1|) and (|5.2[) give rise to an expression for W ac- 
curate through order e 3 . 

The expansion coefficients are determined by inserting 
Eqs. dSTTJ) and {OJ within Eq. (g^J). We begin with W , 
which satisfies 

2a a d a W + (V 2 a + A a a a - 3) W = 0. (5.3) 
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Wc rely on the standard expansions 



0{e 2 ) 



(5.11) 



c a f> — h a a h b b 



ha'b' — -^Ra'c'b'd ,(jC ^ 



(5.4) 



and 
A n = h a ' 



A a , - A a , ]c ,a c ' + -A a , Wd ,<j c ' a d ' + 0(e 3 ) 



(5.5) 

in which a a b '■= D a a b and h a a is the operator of parallel 
transport in the three-dimensional space. From the first 
equation we get 



V 2 a = 3 - li? cM ,a c V + liWi^W + C 



and the second equation gives rise to 



4 ) 

(5.6) 



A a a a = -A d a c + A c ,\ d ,a c a 

2-'c'\d'e' 



1 



A c , ld , e ,a c 'a d 'a e ' +0(e*) (5.7) 



because h a a a a = —o a . Making use of the identity 
cr c a a a = <j c we also find that 

a a d a W a = W°,a c ' + W dd ,a c 'a c ' 

+ lw°, d , e ,a c 'a d 'a e ' +0(e i ). (5.8) 



Making the substitutions within Eq. (|5.3|) and equating 
each expansion coefficient to zero, we eventually arrive 
at 



w°, 

rv a' 


- ^A a , 


v ' a'b' 


= - l -A 
2 


W°,u, , 
v ' a'b' c' 





1 



a'\b' 



(a'\b'c') 



-A„>A h 



;Ra'b', 



(5.9a) 
(5.9b) 



-A( a iA b i\ c >) 



1 



- A a i At,/ A c > 



+ T^(a'-Rft'c') - ^R(a'b'\c')- (5.9c) 

It should be noted that since A a > is the gradient of a 
scalar function, A a ,\ b r — Ar a t\ b ty 

We next turn to W\, which satisfies the differential 
equation 

2a a d a W 1 + (V 2 a + A a a a -l)W 1 = -(V 2 W + A a d a W ) . 

(5.10) 

The left-hand side of the equation is computed with the 
same methods as for the previous computation. For the 
right-hand side we make use of the results a a b , = 5 a b , + 
0(e 2 ), a a ' b = -h\ + 0(e 2 ), and V 2 cr a ' = -\R a ' c ,o c ' + 
0(e 2 ) to obtain 

V 2 W a = h a '»'w% + (-\wlR a ' d + // 6 V a VV 



and 

A a d a W = -A a 'w°, - (A a 'w°, c , - W°,A a ' lc ,)a c ' +0{e 2 ). 

(5.12) 

Making the substitutions within Eq. (|5.10[) . equating 
each expansion coefficient to zero, and simplifying the 
results with Eq. (|5.9[) . we eventually arrive at 



IF 



1 = -A a L + -A a 'A a -R', 

4 |a 8 12 ' 



(5.13a) 



W\ = --A' 



\c' a' 



gA c A d \ a , + -A c \ c ,A a , 



+ ±A c 'A c ,A a , - ^R'A a , + ^R {a ,, (5.13b) 

in which R' stands for the Ricci scalar evaluated at 
x' . The expression for W\, was simplified by invoking 
the contracted Bianchi identity R c a ,\ c , = \R\ a ' as well 
as Ricci's identity to write A a ) c c , + A c ^ a , c , + A c ^ c , a , = 

3A c \ da , + 2R a > c >A c ; recall that A a i\ c i = A d \ a , because 
A a > is the gradient of a scalar function. 

The local expansion of the Green's function is therefore 



G 



H, local 



1 



+ \w% c ,a a 'o b '^ +0(6 4 ) 
+ 2a\w 1 + Wa,a a ' + 0(e 2 ) 
+ 0(a 2 )\, 



(5.14) 



with the expansion coefficients listed in Eqs. (|5.9[) and 



B. Electromagnetic field 

We wish to express the three-dimensional Green's func- 
tion G^(x, x') as a local expansion about the base point 
x' . Once more we rely on the results from the scalar 
case, which we directly import after implementing the 
substitution A a -A a . 

The local expansion of the electromagnetic Green's 
function is 

G"/°™\x, x') = -^={l + W Q a ,a« + \w%a* V 
+ lw°, b , c ,a a 'a b 'a c ' + 0( £ 4 ) 



2(7 



W 1 + Wl,u a ' + 0(e 2 ) 



+ 0(a 2 )\, 



(5.15) 
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with VI. LOCAL EXPANSION IN CONFORMAL 

FORMULATION 

A. Scalar Held 

W°, = -^A a >, (5.16a) The local expansion of Eq. (|5 . 14[) applies to the 

^ ^ i Hadamard representation of the Green's function 

W% = -A a , w + -A a ,A v + -R a ,y, (5.16b) G 3 (x,x' ) defined by Eq. (|3.10[) , but it applies just as well 

to the conformally related Green's function G^(x, x') de- 



2 " |U 4 6 
1 . 3 . 



W a , b , c , = — 2^( a 'l b ' c ') — 7^(a'^&'|c) — g^4a'^6'^4c' hncd by Eq. (|3.13[) ; in this case one simply inserts the 

^ ^ conformally related quantities (such as A a , <r a , and 

_ -^A( a rRy c t} — -R( a 'b'\c>) (5.16c) R a 'b') in place of the original quantities (such as A a , 

a a , and R a 'b')- As we shall now show, the expansions 
are then related by Eq. (|3.17p . 



and 



G^ local 0z, x') = n(x')G^ ocal (x, x'), (6.1) 

a conclusion that guarantees the consistency of the two 
approaches to the local expansion. Thus, a local expan- 
sion formulated in the original space, and a local ex- 
pansion formulated in the conformally related space, will 
W 1 — __ J 4«' ( _|_ —Jf'Aat — R' (5 17a) produce the same Green's function, apart from the fac- 



tor of £l(x') that appears in the relationship between the 



W\ = —A c ' — —A c ' A ii i + —A c ' ,A i Green's functions. 

s ' ca 8 ca 8 ' c This conclusion can be verified by straightforward com- 



—A c ' A i A i + —R' A i + —R (5 17b) putation, making use of the well-known relations between 
16 c a 24 a 24 ' a conformally related quantities. These include 



A a = A a + B ai (6.2a) 

h ab = fl~ 2 h abl (6.2b) 

f to = T f L + 5\B C + S a c B b - h bc B a , (6.2c) 
R\cd = R a bcd + 5 a c D d B b - S a d D c B b - h bc D d B a + h bd D c B a + S a c B b B d - S a d B b B c 

- S a c h bd B m B m + S a d h bc B m B m - h bc B a B d + h bd B a B Cl (6.2d) 
Rab = Rab + D a B b + h ab D m B m + B a B b — h ab B m B m , (6.2e) 

R = n 2 (R + AD m B m - 2B m B m ) , (6.2f) 

D c R ab = D c R ab + D c D a B b + h ab D c D m B m + 2(B a D b B c + B c D a B b + B b D c B a ) 

- (h ac B m D m B b + h bc B m D m B a + 2h ab B m D m B c ) + 2h ab B c D m B m 
+ 2R ab B c + R ac B b + Rb c B a — h ac R bm B m — h bc R am B m + 4B a B b B c 

- (h ac B b + h bc B a + 2h ab B c )B m B m , (6.2g) 
D a R = Q 2 (p a R + 4D a D m B m - AB m D a B m + 8B a D m B m + 2RB a - 4B a B m B m ^) , (6.2h) 

in which B a := d a lnf2, and where all indices on the right-hand side are raised with h ab . They include also 

~d a' i na' a' b' c' i ^ no' a' h' c' d' i nfl' a' b' c' d' p' i AV,-6\ (o n\ 

a =a +-b b , c ,a a a + -S Wdd ,o a a a + — b Vc , d , e ,o a a a a + 0{e ), (6.3) 

an expansion of a a :— h a b D b /a in powers of a a := h a b D b ><j, in which a is half the geodetic separation in the 
conformally related space. The expansion coefficients are given by [3l[ 

S\, c , = 25 a ' b ,B c , - h b , c ,B a ', (6.4a) 
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S", 



b' c' d' 



-26 a b ,D c ,B d , + h b , d D d ,B a + 46 a b ,B c ,B d , - 5%, h c , d , B rn , B m - 2g Vd B a B d ,, 



(6.4b) 



S a b , c , d , e , = 25 a b ,D d D d ,B d - h Wd D d ,D d B a - l25\,B d D d ,B d + Ah b , d B d ,D d B a + A5 a b ,h dd ,B m D e ,B m , 

— h b i c ih d i e iB m D m iB a + 2h b i c iB a D d /B d + 8S a b ,B d B d >B d — A5 a b <h dd iB d B m >B m 

- 4h b , d B d ,B d B a ' + h b , d h d , d B m ,B m ' B a ' - h Vd R a ' d , m , e ,B m \ (6.4c) 

in which the lower indices Ud , or b' d d' ', or b'c'd'e' are understood to be fully symmetrized on the right-hand side of 
the equations. These relations imply 



a(x,x') = n~ 2 (x')cr{x 7 x') 



1 + Pa>° a ' + \Pa> b >o a ' a b ' + l -P a , Vc ,a a 'a b 'cj c ' + 0(e 4 ] 
2 o 



(6.5) 



with 



P a '=B a ,, (6.6a) 

Pa'V = -%D a ,B v + ^B a ,B b , - -h a 'b'B m 'B m ' , (6.6b) 
3 3 6 

Pa'b'c' = -^D a iD b 'B c i — 3B a /D b /B d + -h a 'b'B m D d B m > + 2B a 'B b >B d — -h a i b iB d B m iB m , (6.6c) 

with the same understanding regarding the a'b' or a'b' ' c' indices on the right-hand side. 

I 



To verify that Eq. (|6.1[) holds, we begin with the con- 
formal formulation of Eq. (|5.14[) , in which we make the 
substitutions listed above. Simplifying, and keeping all 
expansions accurate through order e 3 , reveals that in- 
deed, the end result is Eq. (|5 . 14[) formulated in the orig- 
inal space, except for the overall factor of O(x') that 
occurs in Eq. (|6.1[) . Consistency of the local expansions 
is therefore assured. 

It is natural to ask whether the validity of Eq. (|6.ip 
could be established as an exact relation, instead of as 
an approximate local expansion pursued through order 
e 3 . Defining a W(x, x') by the relation 



W(X,X') ;=n( X >)J^l^W(x,x'), 



(6.7) 



the proof would amount to a demonstration that this W 
is suitable to be implicated in a Hadamard construction 
of the conformal Green's function via Gg = W j\[2~a. 

The proof would involve three essential steps. First, 
the function W, as defined here, must be shown to sat- 
isfy the same differential equation as Eq. (|4.3p expressed 
in its conformal formulation; this property follows di- 
rectly from the fact that W = V2&G3, in which the 
two-point function G3 := Q(x')G\l is known to satisfy 
Eq. (|3.13[) . the conformal formulation of Green's equa- 
tion. (The issue at stake is whether this G3, which is 
defined as W / V2&, is a proper Hadamard representation 
of the conformal Green's function.) Second, W must be 
shown to satisfy the boundary condition of Eq. (|4.4[) ; this 
property follows immediately from the coincidence limit 
of Eq. (|6.5p and the fact that W itself satisfies the bound- 
ary condition. Third, W must be shown to be smooth at 
x = x\ by which we mean that the function must be 
G°° when viewed as a function of x with x' fixed; this 



property ensures that W admits an expansion in pow- 
ers of a as displayed in Eq. (|4.5p . which is known to 
be convergent and unique. The expansion being unique, 
smoothness ensures that the Hadamard construction 



(6.8) 



gives rise to a Green's function that satisfies Eq. (|6.ip 
exactly. 

Evidence that W is smooth through order e 3 was pre- 
sented in the context of the local expansion. Because W 
is known to be smooth, smoothness of W to all orders 
relies on the smoothness of <x/<7, which can only be as- 
sured if the series expansion of Eq. (|6.5p can be proved 
to converge. In the absence of such a proof, we shall 
have to give the exact version of Eq. (|6.ip the status of 
a plausible, but unproved, conjecture. 



B. Electromagnetic field 



In Sec. lVI Al we were able to establish that in the case of 
a scalar field, a local expansion of the Hadamard Green's 
function formulated in the original space, and a local ex- 
pansion formulated in the conformally related space, pro- 
duce the same Green's function, apart from the factor of 
fl(x') that appears in Eq. (|3.40p . In addition, we formu- 
lated a conjecture to the effect that the two Hadamard 
forms may be related by 



&£{x,x') =n(x')G^(x,x') 



(6.9) 



as a matter of exact identity. The methods of Sec. IVI Al 
allow us to make the same statements regarding the elec- 
tromagnetic Green's function. The required computa- 
tions are almost identical, and all the relevant equations 
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can be obtained from the scalar case by making the sub- 
stitution A a -> -A a . 



VII. DETWEILER- WHITING CONSTRUCTION 

A. Scalar field 

In this section we construct the three-dimensional ver- 
sion of the Detweiler- Whiting singular Green's function 
for a static scalar field in a static spacetime. By virtue 
of Eq. ()3.16j) . this can be related to the four-dimensional 
version of Eq. ()1.1[) by 

G\{x,z) := f G s 4 (x,z)dr (7.1) 

in which r is proper time for an observer at rest at the 
spatial position z. The integral can be evaluated with 
the techniques described in Sec. 17.2 of Ref. [![, and we 
have that 

G\{x,x) = ^-U(x,x') + -!—U(x,x") 

-±£v(x,z)dr : (7.2) 

in which x' := z(u) is the retarded point on the (static) 
world line, x" := z(v) is the advanced point, r := a a >u a 
is the retarded distance, r a( j v := — a a nu a is the ad- 
vanced distance, and U(x, z) and V(x, z) are the two- 
point functions that appear in the construction of the 
four-dimensional Green's function. 

To calculate G| we follow the methods of Haas and 
Poisson (HP) [32j . wherein the retarded and advanced 
points are related to a middle point x on the world line. 



But while x was chosen arbitrarily in HP, here we specif- 
ically choose x to be simultaneous with x, so that x and 
x have the same time coordinate. This condition implies 
that f := aa{x,x)u a = 0. 

Following HP we define the world-line functions 

cr(r) := cr(x, z(t)j, (7.3a) 
U(t):=U(x,z(t)), (7.3b) 
V(t):=V(x,z(t)), (7.3c) 

in which x is kept fixed. These functions will all be ex- 
pressed as Taylor expansions about r = r, with f defined 
by x, := z(f). We also define 

s 2 :=g & K & a B =2a(x,x), (7.4) 

the squared geodesic distance between x and x. Notice 
that r = <r(u) and r ac i v = —cr(v), in which an overdot 
indicates differentiation with respect to r. We define 

A-:=u-f, A + :=v-f, (7.5) 

with A_ < and A + > 0; these parameters are collec- 
tively denoted A. 

The A parameters arc determined by writing cr(u) = 
or <j(v) = as a Taylor expansion about f: 

= a + a A + -if A 2 + -HA 3 + — A 4 
2 6 24 

+ _L CT (5)A5 + ( e 6 )j (7 . 6) 

in which a and its derivatives are evaluated at r = f. 
This equation is then solved for A. The derivatives of 
ct(t) are given by 



a = \s 2 , (7.7a) 

cr = cr aU a = 0, (7.7b) 
a = (T a pU a u^ + cr & a a 

= -1 - -Ruauo + J^Ruaua;cr + O 'stt" + 0(e 4 ), (7-7c) 

= —^Ruouavu — Ruaaa + C a d Q + 0(e 3 ), (7-7d) 

<^ (4) = OgfapfvP uV + <j & - M (5u & a B v? + u & u r3 a<) + a &fj (Za & a B + 4u a d' 5 ) + a a d & 

= Ruaua -a 2 + a a d & + 0(e 2 ), (7.7c) 

o- (5) = cr a ^ Sg u a u^u^u S u ? + o &m - s (a & vPv?u l + 6u a cAV + 2u & u^a?J + u & u B v?a l ) 

+ a^^cFaPv? + Qu a Ja? + eFvPa? + Vu^aPiP + uF'vPa?) + a^lQePaP + hu & aP) + a a d a 

= -5a a d a +0(e). (7.7f) 
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These results rely on the standard expansion 
(r a p(x, x) = g a p — -R &fi pp<T^ l o v 

+ ^Ram-k^o* + ( 7 - 8 ) 

which can be differentiated with respect to x a to produce 
expansions for and so on. We use the HP notation 
for the components of the Riemann tensor; for example 

Ruaua ■= RafLfiv ua < jllu ^< jV i Ruaaa ■= R & ^U a CJ^ U y , 

and R urT ua:u '■= Raf L fiv:\ ua ' J,lu ^' J ' yu> We have defined 



, (7.9) 



and so on, and used the identities «„a M = 0, tt„a' 1 



—a := a,,a^, and 



-3a, ,a^. 



Substitution of these expansions within Eq. (|7.6[) and 
solving for A returns an expansion of the form 

A = Aie + A 2 e 2 + A 3 e 3 + A 4 e 4 + <3(e 5 ). (7.10) 

The explicit expressions for Ai, A 2 , A3, and A4 are too 
large to be displayed here, but we may mention that 
A+ = s and A^ = — s. 

With A determined, r and r a( j v can be calculated as 
Taylor expansions. Since r = a(u) and r a dv = —&(v), we 
have that 

r = aA_ + -aA 2 + -a^A 3 + —a^A 4 
2 " 6 24 

+ C>(e 5 ), (7.11a) 



?~ach 



= -aA + - iaAi - V>Ai - itr^Al 



0(e 5 



6" ~+ 24 + 

(7.11b) 



At leading order r — s + 0(e 2 ) and r a d v = s + 0(e 2 ), but 
the complete expansions for r _1 and r~^ v are too large 
to be displayed here. 

Expressions for U{x,x') and U{x,x") are obtained in 
a similar way. We write 

U(x, x') = U + U A- + ^UA 2 _ + qUA 3 _ 

+ 0(e 4 ), (7.12a) 
U(x, x") = U + UA+ + ^UA 2 + + ^UA 3 + 



+ 0(e 4 ), 



(7.12b) 



in which U(t) and its derivatives are evaluated at r = r. 
These quantities are given by 

U = 1 + -^R aa - ^iW + 0(e 4 ), (7.13a) 
J7 = U- A u a 

= -xRua + 7^R<?cr-u ~ TTr-Riier:^ + 0(e 3 ), (7.13b) 

24 12 

[/ = U. M u & uP + C/ ;a a a 



6 6 6 12 

+ 0(e 2 ), (7.13c) 

= + + 0(e). (7.13d) 

To evaluate the tail integral we expand V(t) as V(f) + 
(t — f)y(r) + 0(e 2 ) and integrate with respect to r be- 
tween m = f + A_ and v = f + A + . The result is 



y(x, z) dr = F(A+ - A_) + ^V{A\ - A 2 _) 



+ 0(e 3 ), 



(7.14) 



in which V := V(x,x) and V := V ;S u a . These are given 
by the expansions 

V = ±R-±R., a a & + 0(e 2 ), (7.15) 

and 



V = -R. A u & + 0(e). 



(7.16) 



To obtain Eq. 15[) we rely on standard expansion tech- 
niques. The two-point function is required to satisfy the 
wave equation DV — as well as the light-cone equation 

V ]a a a + \{a<* a -2)V = \uU, (7.17) 

which is evaluated at er(x, x) = 0. The solution is ex- 
pressed as an expansion 



V{x,x) = J^V n (x,x)a n , 



(7.18) 



n=0 



and the wave equation gives rise to a sequence of equa- 
tions which determine V n from V n -i \ the light-cone equa- 
tion determines Vq. Because a = 0(e 2 ), V = Vq to order 
e, and this can be obtained by inserting the expansion 



V = V" + Vgo" + Q{e 2 ) 



(7.19) 



within the light-cone equation. We use the fact that 
a a a = 4 + 0(e 2 ), and to compute DU we start with 
Eq. (|7. 13b|) and rely on the expansions 

< = sT a + 0(e 2 ), g%p = O(e); (7.20) 
we eventually arrive at 

OU = \r - \r.. & u & + 0(e 2 ). (7.21) 
6 6 



The end result of the computation is Eq. (|7. 15|) . 

Putting all the ingredients together, we eventually ar- 
rive at the following expansion for Gf : 



13 



with 



3 1 1 



a/3 ~ 4 



6 



l g/3 



(7.22) 

(7.23a) 
(7.23b) 
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a/37 



ipdry — —a a u' J 'u v Rp i p D ^ + —a a Rp^ 



-u ft u u Rp i&j? p. ; y — —R 5 p.^ 



(7.23c) 



and 



^ = -\^a fl + ^v?v?Rp,* - ^R, 



(7.24a) 



5 1 1 1 

— a M a p a Q + -a Q + -a & u^u u R^ - — i?a e 

It) o o ^4 



1 



fa D u X R 



1 



1 



upas t ^ atl Rfj.a + -^2 U ^ U 



24 ^nftOtT 24 



(7.24b) 



The actual expression for V'g^ is obtained from what ap- 
pears above by symmetrizing over all three indices; this 
operation was suppressed to keep the notation unclut- 
tered. 

Noting that the vector a a has a vanishing time com- 
ponent when x and x are simultaneous events, we may 
re-express Eq. (|7.22[) as 

1 



^+V^ s + 0(e 2 ) 



(7.25) 



And with the results derived in Sec. UH the expansion 
coefficients become 



^ = -^815+^-^ + ^. 

^llc = \ A {a\bc) - \ A {a A b\c) + % A a A b A < 
+ -7 A {aRbe) ~ T^(a5|e)- 



and 



y/ = -A\ + -A a A- a - —R, 



(7.26a) 
(7.26b) 

(7.26c) 

(7.27a) 



1 - 1 

— i?A a + —R\a- 



24 



24 



(7.27b) 



Comparing Eq. (|T2"6")) with ([5T5]). and Eq. (j7T2"T)) with 
(|5.13p . we observe that the expansion coefficients of Gf 
and Gg are in precise agreement. This allows us to con- 
clude that 

Gl(x,x)=G%(x,x) + 0(e 3 ) (7.28) 
for a static spacetime. 

B. Electromagnetic field 

We next turn to the three-dimensional version of the 
Detweiler- Whiting singular Green's function for a static 
electromagnetic field in a static spacetime. By virtue of 
Eq. (|3.44p . we have that the vector potential of a point 
charge e situated at z is given by 

${x) = -eN(z)C%(x,z), (7.29) 

with Gf(x, z) denoting the three-dimensional version of 
the Detweiler- Whiting electromagnetic Green's function. 
And according to Sec. 18.2 of Ref. [l|, we have that the 
vector potential is given 

*a(*) = ^CVOz.zV + -^U aP 4x,x") U P" 
2r 2r adv 



V afl (x, z)u^ dr, 



(7.30) 



in which U afjt (x, z), V a)1 {x, z) are the two-point functions 
that appear in the construction of the four-dimensional 
Green's function. 

To calculate $f and obtain Gf we once more follow the 
methods of Haas and Poisson (HP) (32|, as outlined in 
the scalar case. We thus define the world-line functions 



cf(t) := a(x, z(t)), 
U a {r) :=U atl (x,z{T))u^( T ), 
V a (r) :=V aiM (x,z(T))u^( T ), 



(7.31a) 
(7.31b) 
(7.31c) 



in which x is kept fixed. These functions are all scalars 
with respect to their dependence upon z(t). As in the 
scalar case they are expressed as Taylor expansions about 
r = f, at which z = x, and the results are converted into 
explicit expressions for r, r a( jv, U a fjiu^ , U a p"U^ , and 
the tail integral. The results for r and r a( j v appear in 

Eq. frm 

To compute U a piu° and U a /3"U^ we write 



{7 Q A_ + l-U a A 2 _ 



1_ 



+ 0(e% 



(7.32a) 



UapnuP" = U a + C/ Q A^ 
+ 0(e% 



^U a A 2 + + \tJ a A s + 



(7.32b) 
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in which U a (r) and its derivatives are evaluated at r = r. 
These quantities are given by 



14 



1 + ^Raa - ^Raa;a + 0(e 4 ) ) 



a:a;/3 
Qaa. 



I* 



i 



i 



24 



12 



1 



12 



(7.33a) 



(7.33b) 



U a = U a& .- M u & u p v? + U a& . B (u a a p + 2a a u") + U a& a & 



q R uucr-.u r,R uacr R aua u ( R Ruu + qRu<j;u ir) -Rtm;cr + a R 



12 



+ a a ( ^R ua )+a & + 0(e 2 ) 



(7.33c) 



U a = U a& .fogU S uPu^u* + U aa . f ^(3a a u^iP + 2u & a^v? + u s u P a^) + I ,,„ : A : '>"" " ' + 3o a u^ + u 5 ^) + U aBt a & 



-R a uau + u a -R au + -R uu , u )+a a [ -R uu + a" + 0(e) 



(7.33d) 



r 



The expansions involve components of the Riemann ten- which leads to 
sor such as R a uua := R a -p-u fl u l3 <j L ' and components of 

the Ricci tensor such as R aa := R a sa a a^ . They involve V a = g aa 

also g a a (x,x), the parallel propagator from x to x. To 
arrive at these results we rely on the expansion of the _j_ ^, 

two-point function U a a (x,x) given by 



U% = g\ (l + y 2 R-^°* - ^R^ot'o* + 0(e 3 ) 

(7.34) 

this is derived, for example, in Appendix B of Ref. [33| ■ 
Another useful expansion is 



g a _ - = g% ( -R^_-a^ - -iT_ 5 - -^v" + 0(e 3 ) 



}_ pa , }_ r>a j \_ n a !A" 

2 f 4 ' 19 



^-21^ 



_ }_ r>a j _ pa _ _ no 



+ 2 ^« a ^ + ^« a J R + 0(e) 



(7.38a) 



(7.38b) 



Here we make use of the notation R a a := i? Q ^a' 3 , R a u . a := 

R a g.-V>t ', and i?"^'^ ■= R a ^-- ,tl u l3 u 1 . In addition, 

is the Ricci scalar evaluated at x, and -R ;u := R.avF . 
To obtain Eq. (|7.37p we rely on standard expansion 



(7.35) 

These are differentiated repeatedly with respect to x a , 
and the results are inserted within the expressions for techniques. The two-point function is required to satisfy 
U a and its derivatives. the wave equation 

To evaluate the tail integral we expand V a (r) as 



Vq(t) + (t — r)V r Q (r) + 0(e 2 ) and integrate with respect 
to t between u = f + A_ and v = f + A+. The result is 

f 7„ M «"dT = V a (A+ - A_) + iy«(A^ - A 2 _) 

+ 0(e 3 ), (7.36) 



□n - = o 



as well as the light-cone equation 



(7.39) 



+ - 2)V% - \{UU a a - R a p U^), (7.40) 



in which V a = V a& u & and V a = V aa .- p u & u B + V" QQ a a . To pressed as the expansion 
compute these quantities we rely on the expansion 



V a = q a 



Ry_ - tp.R 
i „ 



which is evaluated at <r(x, x) = 0. The solution is ex- 
expansion 

V%(x,x) = J^V£ a (x,x)a n , (7.41) 



12 a M/3 



n=0 



and the wave equation gives rise to a sequence of equa- 
(7.37) tions which determine V" a (x,x) from V"_ la ; the light- 
cone equation determines Vq & . Because a = 0(e 2 ), 
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V % = V a a to order e, and this can be obtained by in- 
serting the expansion 



V Q a a =9% A\ + A\-^ + 0{e 2 ) 



(7.42) 



within the light-cone equation. We use the fact that 
a p — 4 + 0(e 2 ) and to compute OU a a we rely on the 
expansion 



Oi OL P 

9 &-ft = 9 -,9 a 



we eventually arrive at 



\5\R-\5\R,^ + V, w „/V + 0(e 2 ) 



(7.44) 

The end result is Eq. (|7T37)) . 

Putting all the ingredients together, we eventually ar- 
rive at the following expansion for $>^(x): 



+ ^W S ^^ + 0(e 4 ) 



a a + 0(e 2 ) 



(7.45) 



(7.43) with 



and 



<t>l 



"A J 




1 




2 A 
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— a 

8 



3 - - 1 



-a & R^ -R & p-^ 



(7.46a) 
(7.46b) 

(7.46c) 
(7.46d) 



( - + -^u»u v R^ -^Rj+±a- x + ^R- x -^, 

- ^-cfiajiCbo, + ia g + \a a u^u v Rj^v - -^Ra & + \u^a" u p R^p & + —aPRjxa 

10 o o Z4 o Iz 



12 
1 



24 
1 



6 



12 



2 ' 6 
1 



(7.47a) 



A/Ipa 1" ^ utl R\fi a a — R\p,-a- 



(7.47b) 



The actual expression for <^j| a ^ is obtained from what 
appears above by symmetrizing over the last three in- 
dices; this operation was suppressed to keep the notation 
uncluttered. 

From Eq. (|7.45j) we wish to obtain a more explicit ex- 
pression for $j , and this requires a computation of the 
operator of parallel transport. Our considerations near 
Eq. (|2.5[) imply that its components are given by 



N(x) 



a ia 

9b = h k> 



{7.4 



in which h a r is the operator of parallel transport in the 
three-dimensional space; the mixed components c/' a and 
g a r vanish. Noting that the vector a a has a vanishing 



time component when x and x are simultaneous events, 
and making use of the results derived in Sec. [TT1 we may 
re-express Eq. (|7.45[) as 



= --N{x)\ 1 + 4>la a + 
1 



with 



6 



4> l + 4>W a + 0(e 2 ) 



:A- n 



(7.49) 



(7.50a) 
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4^315 +4^^ + ^85' ( 7 ' 50b ) 



3 1 

• J abc ~ 2^(5|bE) ~ 1 A {a A b\c) + gAaA h A £ 



-A(aRbc)- J R (ab\c)- (7.50c) 



and 



(7.51a) 



o^ C |ea — -A c A £ \ a — -A c \ 5 A a + —A c A E A a 



\ca g 
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24 
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(7.51b) 



From Eq. ([739)) and (f7T49|) we see that the 
three-dimensional Green's function involves the ratio 
N(x)/N(x). This can be expressed as an expansion 
about x = x by making use of the generalized Taylor 
series 

N(x) = N(x) - N la a a + \N^a a a l - ^A^aVV 



+ 0(e% 

which leads to 
N(x 



(7.52) 



N(x) 



= l-A u a a + -(A- afb + A- a A h )a a a b 
- \{A- a \bc + 3A- a A- b{5 + A- a A- b A e )a a a b a s 



+ 0{e i ). 
With this we finally arrive at 



Gl{x,x) = \{l + ^ + \^a l 

1 _ 

/,0 _a b _c 



(7.53) 



with 



V4 - \A- a \i + \A- a A b + ii? aS , 

€fc = -\ A {a\bc) ~ \ A {a A b\c) - \ A -aA b A e 
~ J A (aRbc) - J R (ab\c)- 



(7.54) 

(7.55a) 
(7.55b) 

(7.55c) 



and 



(7.56a) 



1 



1 



1 



i'a — % AC \ca -^A C A e \ & + % A °\c A a ~ jQA C A c A a 



24 RAa + 24^ 



(7.56b) 



We notice that the expansion coefficients can be obtained 
from Eqs. (|7.26j) and (|7.27j) by making the replacement 
A a —> —A a ; this was expected since Eq. (|3.35|) for the 
electromagnetic Green's function differs from Eq. (|3.10[) 
for the scalar Green's function by the sign of A a . 

A comparison between Eq. (|7.54[) and Eq. (|5 . 1 5[) allows 
us to conclude that 



Gl(x,x) = G^(x,x) + 0(e 3 ) 
for a static charge in a static spacetime. 

VIII. EQUALITY OF G| AND G": A 
CONJECTURE 

A. Scalar field 



(7.57) 



The result of Eq. (|7.28[) suggests that the equality 
between the Hadamard and singular Green's functions 
might be exact, holding to all orders in e. We re-express 
Eq. OJ as 



Gf(sc,sc) = -W s (x,x) 
s 



(8.1) 



with 



W s (x,x) 



-U(x,x') + -U(x,x") 
r r 



— s V(x, z(r)) dr 



• 2) 



and conjecture that W s = W H , where W H is the two- 
point function introduced in Eq. (|4.2[) . We recall that 



,s 2 := 2a(x,x) is the squared geodesic distance between 
x and the simultaneous event x, x' := z(u) is the re- 
tarded point on the (static) world line, x" := z(v) is 
the advanced point, r := a a 'U a is the retarded distance, 
^adv := —<J a "U a is the advanced distance, and U(x, z), 
V(x, z) are the two-point functions that appear in the 
construction of the four-dimensional Green's function. 

As in Sec. IVI Al above, a proof of equality would in- 
volve three essential steps. First, the function W s must 
be shown to satisfy the same differential equation as W H , 
as displayed in Eq. (|4.3[) ; this property follows immedi- 
ately from the fact that G§ is known to satisfy Eq. (|3.10D , 
just like G3. Second, W s must be shown to satisfy the 
boundary condition of Eq. (]4.4[) : this property was estab- 
lished previously and can be seen directly from Eq. (|7.25p . 
Third, the function W s must be shown to be smooth at 
x = x, by which we mean that the function must be C°° 
when viewed as a function of x with x fixed; this property 
ensures that W 5 admits an expansion in powers of a as 
displayed in Eq. (|4.5p . which is known to be convergent 
and unique. The expansion being unique, smoothness 
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therefore ensures that W s = W H . Evidence that W s is 
smooth to order e 4 is provided by Eq. (|7.25[) . 

Of the ingredients involved in the make-up of W s , the 
two-point functions U(x, z) and V(x, z) are known to be 
smooth, but s, r, r a d v , u, and v are not. Nevertheless, 
we conjecture that the combinations 



s/r, s/r adv 



s(v — u) 



(8.3) 



are in fact smooth at x — x. The first two are directly 
involved in Eq. (|8.2[) , and the third one also is involved 
by virtue of the mean-value theorem, which allows us to 
write the integral as 



V(x, z) dr = V(x, x*)s(v — u), 



(8.4) 



with x* := z(t*) (u < t* < v) representing a middle 
point on the world line. Establishing that s/r, s/r a( j v , 
and s(v — u) are smooth is sufficient to prove that W s 
itself is smooth. 

Some insight can be gained by examining these quan- 
tities in Fermi normal coordinates (t,x a ) attached to the 
static world line. With results collected from Sec. 11 of 
Ref. [l[, we have that 



\J S a bX a X b 



1 + ^a a x a 



-(a n x a ) 2 



-a t s 



1 



R tatb x a x b + 0(s 3 ) 



6 

r + <3(s 4 



^*adv 

U = t — S 
1 



1 3 

1 - -^a a x a + -(a a x a 



24 



^RtatbX a X b + 0(S 3 ) 



t + s 

1 

~ 6 



1 3 

-a a x a + - (a a x a 

2 8 V 



24 



a t s 



Rtat b x a x b + 0{s 3 ) 



3.5a) 

5.5b) 
3.5c) 

5.5d) 
3.5e) 



in which a a , at, and Rtatb respectively represent compo- 
nents of the acceleration vector, its covariant derivative, 
and the Riemann tensor evaluated on the static world 
line, at which x a = 0; terms involving a a were discarded 
because these components vanish for a static world line 
in a static spacetime. These results reveal that s,r, r a d v , 
u, and v arc indeed not smooth at x a = 0. But they do 
show that r/s, r^/s, and 



s(v — u) = 2s* 



1 - \a a x a + ^(a a a; a ) 2 



1 ■ 2 
— ats 
24 



-RtatbX a X b + 0(S 3 ) 

6 



(8.6) 



are smooth to leading orders in an expansion in powers 
of x a . 



We now proceed with a sketch of what might consti- 
tute a general proof. The method of proof relies on formal 
power series, which are all assumed to converge in a suf- 
ficiently small domain. This rather strong assumption is 
the main limitation of our argument, and the reason why 
we present it as a conjecture and not a proof. It would 
be desirable to either establish the convergence property, 
or to devise an alternative method of proof. This shall 
be left for future work. 

We return to Eq. (|7.6j) and observe that the odd terms 
in the expansion vanish by time-reversal invariance: a, 
a, and all other odd derivatives of <t(t) must vanish 
on a static world line in a static spacetime. We recall 
that ct(t) := o~(x, z(t)) with x fixed, and state that each 
derivative of <t(t) is smooth at x = x. Equation (|7.6j) 
can therefore be written as 



,s 2 = A 2 £p„(A 2 y 



n=0 

2 



f(- 



r (2n+2) 



(8.7a) 



.7b) 



(2n + 2)! 

in which a bracketed number attached to a indicates the 
number of differentiations with respect to r; each expan- 
sion coefficient p„ is smooth at x = x. Time-reversal 
invariance implies that A± = ±V A 2 := ±A, and the 
expansions of Eq. (|7.11|) can be expressed as 



r adv = A^ g „(A 2 )", 



1 



r (2n+2)^ 



with q n smooth at x 
have that 

r r adv 



(8.8b) 

(2?i + l)! v v ; 

x. Combining these results, we 



2\n 



£~=pgn(A 



(8.9) 



As stated previously, each sum in this expression is as- 
sumed to converge for A 2 sufficiently small. 

We now wish to reverse the expansion of Eq. (|8.7|) . 
According to Sec. 3.6.25 of Ref. |3J], if y — ax + bx 2 + 
cx 3 + • • • , then x = Ay + By 2 + Cy 3 + ■ ■ ■ with aA = 1, 
a 3 B = —b, a°C = 2b 2 — ac, and an algorithm is known to 
generate all remaining expansion coefficients. The power 
series can thus be reversed when a =/= 0. In our case 
a = pq = —6 is indeed nonzero, and a -1 is smooth at 
x = x. The reversed series can then be written as 



n=0 



(8.10) 



for some coefficients a n that are known to be smooth at 
x = x. Because s 2 is itself smooth, the assumed conver- 
gence of the sum for sufficiently small s 2 ensures that A 2 
is smooth at x = x. Making the substitution in Eq. ()8.9[) . 
we find that r/s and r a dv/s can be expressed as 



£~oM* 2 ) n 



(8.11) 
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for some coefficients b n and c„ that are smooth at x = x. 
This reveals that r/s and r ac i v /s are smooth at x = x. 
We next turn to s(v — u) = 2sA, which is given by 



s(v - 



5.12) 



and is also seen to be smooth at x = x. 

With the stated assumption on the convergence of for- 
mal power series, we have shown that r/s, r a d v /s, and 
s(v — u) are all smooth at x = x. This implies that W s 
is smooth, and establishes the statement that Gf and G3 
are strictly equal. 



B. Electromagnetic field 



The simplicity extends to the two-point function 
U(x,x') that enters the Detweiler- Whiting construction. 
We may show, in particular, that U has no dependence 
on the time coordinates, so that 



U = U(x,x'). 



(9.3) 



This statement is a consequence of the defining properties 
of the two-point function (see Sec. 14.2 of Ref. [lj), that 
it must satisfy the differential equation 



2a a d a U + (a a a -4)U = 



(9.4) 



in the ultrastatic spacetime, together with the coinci- 
dence limit U(x',x') = 1. With the stated properties of 
the world function, this becomes 



The result of Eq. (|7.57l) suggests that the equality 
between the Hadamard and singular Green's functions 
might also be exact in the case of the electromagnetic 
field. A proof of this statement would involve the same 
steps as in the scalar case, and the modifications required 
for the electromagnetic Green's functions are too modest 
to merit a separate discussion. As in the scalar case, the 
essential element is the proof s/r, s/r a( jv, and s(v — u) 
arc all smooth at x = x. If this can be established, then 
we can claim immediately that G| and G3 are indeed 
equal to all orders. 



IX. EQUALITY OF G| AND G3 FOR 
ULTRASTATIC SPACETIMES 

A. Scalar field 

In this section we return to the theme explored in 
Sec. IVIIII and provide a complete proof of equality be- 
tween the Hadamard construction G3 and the three- 
dimensional version of the Detweiler- Whiting Green's 
function Gf in the case of ultrastatic spacetimcs. These 
spacetimes have the property that N — 1, so that their 
metric is 



ds 2 



-dt 2 + h ab dx a dx b 



(9.1) 



a special case of Eq. (|2.1j) . The geometrical quantities 
associated with ultrastatic spacetimes can be obtained 
from the equations displayed in Sec. HI] by setting A a :~ 
d a lnN = 0. 

The geodesies of ultrastatic spacetimes are described 
by the equation t(A) = £(0)+t(0)A, in which A is an affinc 
parameter and an overdot indicates differentiation with 
respect to A, as well as the statement that x a (X) describes 
geodesies of the spatial metric h a b. This implies that the 
world function is necessarily given by 



<t{x, x') = --{t - t 



l\2 



a 3 (x,x'), 



(9.2) 



in which a 3 (x, x') is the three-dimensional version of the 
world function, defined with respect to the spatial metric. 



(t - t')d t U + oldJJ + l - (a 3 a a -3)U = 0. (9.5) 

The differential equation can be integrated along any 
spacetime geodesic that originates at x' . We may, in par- 
ticular, choose a time-directed geodesic with no spatial 
displacement, such that t(X) = t' + A and x a (X) = x a . 
For such a geodesic we have that 0-3 = and a 3 a a = 
a 3 a a (x' , x') = 3, and the differential equation reduces to 
(t — t')dtU = 0. This implies that the two-point func- 
tion cannot depend on t, and since its dependence on t' 
can only be through the combination t — if, it cannot 
depend on t' . We have, therefore, established the stated 
property. 

The absence of a dependence upon t implies that the 
two-point function satisfies the purely spatial differential 
equation 



2^3 a 



Z)U = 



(9.6) 



together with the boundary condition U(x',x') 
1 . These are precisely the defining relations for the 
Hadamard function Wq(x, x'), as stated in Eqs. (|4.6j) and 
(|4.7p . We conclude, therefore, that 



U(x,x') = W (x,x') 



■7) 



for ultrastatic spacetimes. 

Next we turn our attention to the two-point function 
V(x,x'), and prove that it admits the expansion 



00 

V(x,x') = Y,Vn(x,x')a r ' 

n=0 



(9.8) 



in which the coefficients V n are smooth and time- 
independent; the expansion involves the four-dimensional 
world function, and it is known to converge within a suffi- 
ciently small neighborhood of x' . The proof of the state- 
ment relics on the recurrence relations satisfied by the 
expansion coefficients [30| . 



a a d a V + ^{a a a -2)V Q = laU 



(9.9) 
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when n = 0, and 
<r a d a V n + ±{* c 



2n-2)V n = -—DV„-i (9.10) 



when n > 0. 

We begin with an examination of Vq. In ultrastatic 
spacetimes its differential equation becomes 

(f - t')d t V + <T$d a V + ±{a 3 a a -1)V = \v 2 U. (9.11) 

Once more this equation can be integrated along any 
spacetime geodesic that originates at x' , and once more 
we choose a time-directed geodesic. In this case we have 



(t - t')d t V + V = iv 2 C/ 



^R(x'), (9.12) 



in which R(x') is the spatial Ricci scalar evaluated at x', 
obtained from the known expression for V 2 C/ evaluated 
in the coincidence limit (see Sec. 14.2 of Ref. [l[). The 
general solution to this equation is Vq = j^R(x') + c(t — 
i') -1 where c is a constant, and we see that Vq fails to 
be smooth at x = x' unless c = 0. We conclude that Vq 
cannot depend on time. 

Turning next to V n , we proceed by induction. We as- 
sume that V n -\ is known to be time- independent, and 
prove that V n must in turn be time- independent. We 
begin with the differential equation 

(t-t')d t V n +(T%d a Vn + ^(<T 3 a a + 2ri-l)Vn = "^VVn-l, 

(9.13) 

which we integrate along a time-directed geodesic. The 
equation reduces to 



(i - t')d t V n + (n + l)V n = -t^V 2 K-i 

In 



(9.14) 



and we find that the general solution contains a term 
c(t — t')~( Tl+1 ) that fails to be smooth at x = x' unless 
c = 0. This allows us to conclude that V n cannot depend 
on time, and we have established Eq. (|9.8|) . 

We may now demonstrate the equality of the Green's 
functions. The four-dimensional version of the Detweiler- 
Whiting singular Green's function is 



1 



1, 



with At = t — t' , we find that the integral becomes 
U{x,x') 1 



G\{x,x') = 



V{x,x')dAt. (9.18) 



In this we insert Eq. (|9.8[) . integrate term by term using 



a n dAt = I — / (At 2 - 2a 3 ) n dAt 
V^(n+ 1 



2 n T(n+ I) 



(9.19) 



and simplify. Our final expression for the singular 
Green's function is 



G s 3 (x,x') 



W s {x lX ') 
x/2^7 



(9.20) 



with 

W 5 {x,x') = U(x,x') 

-Y,jt^h v -^ x ' ){2 ^ )n - (9 - 21) 

n=l y '" 

These equations reveal that G% does admit a three- 
dimensional Hadamard form, and that we may make the 
identifications 



Wq(x,x') := U(x,x') 
as in Eq. (|9~7|) . and 

Wl{x,x') := - ^-%\ Vn - lix ' X,) - 



(9.22) 



(9.23) 



Gl(x, x') = -U(x, x')6(cr) - -V(x, x')Q(<j), (9.15) 



These coefficients satisfy the recursion relation of 
Eq. (|4.6p . as can be seen by invoking Eqs. (|9.9p and 
(|9.10p , and are therefore the same coefficients that appear 
in Eq. (|4.5p . The proof of equality between G\(x,x') 
and G^(x,x') in ultrastatic spacetimes is complete, and 
the calculations have revealed the relationship between 
U and Wq, and between V n and W n . 



B. Electromagnetic field 



in which O is the Heaviside step function and S the 
Dirac distribution. According to Eq. p,16[) . the three- 
dimensional version is 



when N(x') 
factorized as 



Gl(x,x')= I Gl(x,x')dt' (9.16) 
= 1. With U independent of time and a 



i(A< 



IW, ){At 







(9.17) 



The proof of equality between the Hadamard con- 
struction Gg and the three-dimensional version of the 
Detweiler- Whiting Green's function G\ for ultrastatic 
spacetimes proceeds along the same lines as in the scalar 
case. In fact, the calculational details are strictly iden- 
tical, because the two-point functions £// (x^x 1 ) and 
V t * (x,x r ) that are involved in the relevant component 
of the electromagnetic Green's function, 

G s t\x,x') = if//(x,x')<5(a)-^ t t '( a ;,x')0( ( T), (9.24) 
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are strictly identical to their scalar counterparts: [// = 
U and V t l =V. The first equality follows from the gen- 
eral relation Uf = g/lJ — see Eqs. (14.8) and (15.9) 
in Ref. [l| — together with the property that gl = 1 for 
ultrastatic spacetimcs. The second equality follows from 
the fact that if V£ is expanded as 



a / j n a 



(9.25) 



then the recursion relations satisfied by V n t * are strictly 
identical to those satisfied by V n . The results of 
Sec. IIX A[ therefore, allow us to state that for ultrastatic 
spacetimcs, Gij 1 = Gf in the electromagnetic case also. 



X. SELF-FORCE IN SPHERICAL SPACETIMES 

A. Scalar field 

We consider the self-force acting on a static scalar 
charge q in a static and spherically-symmetric spacctime. 
The metric is written as 



in which a prime indicates differentiation with respect to 
r. The modes with m ^ necessarily vanish. 

The self-force acting on the scalar charge is given by 
F a = q(g al3 + u q m' 3 )V / 3$ r , in which $*:=$- $ s is 
the difference between the actual potential $ and the 
Detweiler- Whiting singular field 3> s ; the regular potential 
is known to be smooth at x = z. In a static situation 
the self- force has a vanishing time component, and its 
spatial components are given by F a = qh ab db& R - In a 
spherically-symmetric spacetime the angular components 
vanish, and the radial component is 



F r = qfod r <P R (r , 



?0,<P0 



(10.6) 



Recalling the spherical-harmonic decomposition of the 
potential, we may express this as 



in which 



(d r <P) e := V *i m (r)y&, 



(10.7) 



(10.8) 



ds 2 = -e 2 ^dt 2 + r l dr 2 + r 2 



; 6»# 2 ), (10.1) 



in which ip and / are functions of r. In this notation 
N = e^ and A r = tp' is the only nonvanishing component 
of the vector A a . The potential $ generated by the point 
scalar charge is a solution to 



V 2 $ + A a d a $ = -ATrqS 3 (x, z), 



(10.2) 



which is obtained from Eqs. (|3.2|) and (|3.21|) . 

In order to integrate the field equation we decompose 
the potential and source in spherical harmonics: 



§{t,6,<t>) = Y,®im{r)Y l 



I'm { 



(10.3) 



im 



and 



5 3 (x, z) = ^6(r - r ) £ Y e * m (9 , cf> )Y im (e, 0), 

r im 

(10.4) 

in which (ro,0o,4>o) represent the spherical coordinates 
of the particle's position z, and f :— f(r a ). With- 
out loss of generality we may place the particle along 
the polar axis (#o = 0) and exploit the property 
Yem(0,<P) = \J (2^ + l)/(47r)<5 mj o of spherical-harmonic 
functions. Substitution within the field equation then 
produces 



1/2 



2 + jj- + rip' 



£(£ + !) 
f 



<I>, 



= —Anq 



^V /2 *(r-ro), (10.5) 



are the multipole coefficients of d r $>, while (9 r $ s )^ are 
those of the singular potential $ s . Recalling the relation 
of Eq. (|3.22[) between the potential and the scalar Green's 
function, we may write this in the form 



F r = q 2 f lim J2 [^T 1 ( d r$) e - {d r G s 3 ) 



(10.9) 



in which G|(a;, z) is the three-dimensional version of the 
Detweiler- Whiting singular Green's function introduced 

in Sec. IVnAl 

The limit in Eq. (|10.9|) can be taken by setting r = 
ro + A, 8 = 8q, <f> = <f)Q, and letting A — > (from either 
direction). With this choice, we shall show below that 



(d r G 



3)i 



A 



B 



C 



D 



+ o(r 3 ), 



(10.10) 



in which the regularization parameters A, B, C, and 
D depend on ro but are independent of £; explicit ex- 
pressions will be presented below. Inserting Eq. ([lO.lOp 
within Eq. (|10.9[) provides a practical method of com- 
puting the self-force by means of a regularized mode sum 
that converges to the correct answer. With the particle 
placed on the polar axis (#o = 0), the multipole coeffi- 
cients reduce to 



2£+ 1 

47T 



$So(r + A). 



(10.11) 



To establish the relation of Eq. (|10.10[) and calculate 
the regularization parameters we follow the method de- 
scribed in Sec. V of Haas and Poisson (HP) [32|, which 
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we adapt to the situation at hand. In HP the motion of 
the particle was geodesic and the spacetime was that of 
a Schwarzschild black hole; here the particle is kept in 
place in any static, spherically-symmetric spacetime. In 
HP the motion was taking place in the equatorial plane, 
and a transformation of the angular coordinates was im- 
plemented to put the particle momentarily on the polar 
axis; here the particle is kept on the axis at all times, and 
the transformation is not required. Following Sec. Ill of 
HP, the singular Green's function of Eq. (|7.25|) is ex- 
pressed as an expansion in powers of the coordinate dis- 
placements w a := x a — x a , in which x := z denotes 
the particle's position. As in Sec. V of HP we express 
the angular separations w e and in terms of functions 
Q := \/l — cos 6, sin</>, and cos</> that are globally well- 
defined on the sphere. In this case of static motion, the 
squared-distancc function introduced in Eq. (5.22) of HP 
reduces to 

p 2 = f^A 2 + 2r 2 Q = 2rt(5 2 + l-cos6), (10.12) 
where 

A 2 



5 2 



2r 2 f 



(10.13) 



with A := w r = r — tq. Following the steps outlined in 
Sec. VE of HP, we obtain an expansion for d r G\ that 
takes the schematic form of 

d r G\ = {d r Gl)_ 2 + (5 r G|)_ 1 + {d r G s 3 ) Q + {d r G s 3 ) 1 

+ 0(e 2 ), (10.14) 

in which a subscript attached to enclosing brackets indi- 
cates the scaling with powers of e. The various terms are 
schematically given by 



{d r G s 3 )_ 2 = M^ 2 (A/p 3 ), 
(d r G s 3 )_ 1 = M_ 1 (l/p) + 0(A 2 /p 3 ) 
+ 0(A 4 / P 5 ), 



(10.15a) 
(10.15b) 



(d r G s 3 ) Q = 0(A/p) + 0(A 3 /P 3 ) + 0(A 5 /P 5 ) 

+ 0{A 7 /p 7 ), (10.15c) 

{d r G s 3 ) 1 = M lP + 0(A 2 / P ) + 0(A 4 /P 3 ) + 0(A 6 /P 5 ) 
+ O{A 8 /p 7 ) + O{A 10 /p 9 ). (10.15d) 

The terms involving the coefficients M_2, M_i, and Mi 
are those giving rise to the regularization parameters; all 
other terms are unimportant. 

The multipolc decomposition of d r G 3 is next carried 
out with the help of Eq. (A19) of Haas and Poisson; be- 
cause the expressions are all ^-independent (by virtue of 
the axial symmetry of the problem), there is no need to 
perform the (^-average described by Eq. (A13). We make 
use of the relations 

,1/2 

(A/p 3 ) e = {l+\)i\- sign(A) + O(A), (10.16a) 



(1/ P )e = - + 0(A), 

TO 

(P)t 



0(A) 



(10.16b) 
(10.16c) 



and arrive at the expression of Eq. (|10.10p with A = 

A/_ 2 / 1/2 r - 2 sign(A), B = Af_i/r , C = 0, and D = 
— MiTq. The detailed computation reveals that 



A 



D 



2? 



/-^signCA), 
3(l + rV), 



(10.17a) 
(10.17b) 



C = 0, 



D 



16r 2 



(10.17c) 

(1 + n/>') - (1 + rip' + 3rV 2 + r 3 ^' 3 



6r 2 ?/;" - 2r 3 ^"')/ + (l + 4r^' + 3r 2 ^")r/' 



+ (l + r^')r 2 /" 



(10.17d) 



in which all functions are to be evaluated at r = ro- 
These are the regularization parameters for a static scalar 
charge in any static, spherically-symmetric spacetime. 



B. Electromagnetic field 

We next consider the self-force acting on a static elec- 
tric charge e in a static and spherically-symmetric space- 
time with the metric of Eq. (|10.1[) . The vector potential 
$t generated by the point charge is a solution to 



V 2 $ t - A a d a <^ t = 4neN(z)6 3 (x, z), 



(10.18) 



which is obtained from Eqs. (|3.29[) and (|3.43l) . 

As in the scalar case we decompose the potential and 
source in spherical harmonics and place the particle along 
the polar axis (9o = 0). Substitution within the field 
equation then produces 



47re 



] j^±le^f~ 1/2 5(r-r ) 1 (10.19) 



in which /o := /(ro), ipo := - 0( r o) J and a prime indicates 
differentiation with respect to r. The modes with m ^ 
necessarily vanish. 

The self-force acting on the scalar charge is given by 
F a = eF R a p u , in which F R a p := F a p - F s % is the 
difference between the actual electromagnetic field and 
the Detweiler- Whiting singular field; the regular field is 
known to be smooth at x = z. In a static situation the 
self- force has a vanishing time component, and in spher- 
ical symmetry its radial component is 



F' 



ee 



-*7o3r$ t R (ro,0o,<fo). 



(10.20) 
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We express this as 

F r = ee -l>°f ^ [{d r $t 
i 

in which 



(10.21) 



21+1 
An 



& U0 (r +A) 



(10.22) 



are the multipole coefficients of d r $t, while (d r ^)e are 
those of the singular potential. Recalling the relation of 
Eq. (|3.44j) between the potential and the scalar Green's 
function, we may write this in the form 

F r = eVo hm ^[e-^-^^^^+^G^J, (10.23) 



in which G\{x,z) is the three-dimensional version of the 
Detweiler- Whiting singular Green's function introduced 
in Sec. IVlTBl 

The limit in Eq. (|10.9[) is taken by setting r = ro + 
A, 6 = O = 0, and <j> = (b = and letting A -> 
(from either direction). With this choice, the multipole 
coefficients of the singular Green's function take the same 
form as in Eq. (110. 10j) . In this case, however, because of 



the different sign in front of A a in the Poisson equation 
for $ t , the regularization parameters are given by 



A = -i/- 1 / 2 sign(A), 



(7 = 0, 



(10.24a) 

(10.24b) 
(10.24c) 



D = * (1 - rib') - (1 - rib' + 3rV 2 - r 3 ib' 3 

16r z l ' v 

+ 6rV + 2rV")/ + (l - W - 3r 2 yL'")rf 
+ (l-ri/'')r 2 H, (10.24d) 

in which all functions are to be evaluated at r = ro. 
These are the regularization parameters for a static elec- 
tric charge in any static, spherically-symmetric space- 
time. The computations that lead to Eq. (|10.24[) involve 
the same steps as those described in Sec. IX Al 
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